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Space-time correlation functions constitute a useful instrument from the research toolkit of 
continuous-media and many-body physics. We adopt here this concept for single-particle random 
walks and demonstrate that the proposed space-time velocity auto-correlation function reveals cor- 
relations which extend in time much longer than estimated with the commonly employed temporal 
correlation functions. We uncover a generic feature of random-walk processes which is the presence 
of a region of strong velocity anti-correlations around the starting point and relate these correlations 
to the density of walkers. To highlight the practical relevance of our findings we outline a setup 
to measure the space-time velocity correlation function of an ensemble of cold atoms evolving in a 
laser-created optical potential. 
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Introduction. Finitcness of velocity is a fundamental 
property of any physical process taking place in space 
and time. This concept was incorporated in the frame- 
work of random- walk theory [l[ ; one of the cornerstones 
of statistical physics that is particularly successful in de- 
scribing diffusion phenomena [11 . If the speed of a walk- 
ing particle is constant, the coupling between the distance 
traveled and the time it takes leads to the confinement 
of the process to a casual cone. Within the cone the 
density of particles is described by the phenomcnologi- 
cal diffusion equation Q. The coupling also regularizes 
fast super-diffusion by removing possibly unphysical di- 
vergences from the momenta of the corresponding pro- 
cesses 0, 

The velocity of a random-walk process can be treated 
as an additional dynamical variable whose evolution is 
itself a random process. The Green-Kubo relation @ 
highlights the importance of the corresponding tempo- 
ral auto-correlation function by connecting its integral 
with the diffusion coefficient of the transport process 
in the original space It is evident, however, that 
for random walks the velocity auto-correlations cannot 
last longer then the time between two consecutive re- 
orientation events of a random walker. 

Following this premise, we next ask whether more ex- 
tended correlations can be detected by unfolding the ve- 
locity correlation function into the spatial domain. Here, 
we answer this question positively by introducing a new 
characteristic that reveals hitherto unnoticed properties 
of random-walk processes. We define the space-time [s-t) 
velocity auto-correlation function for single-particle ran- 
dom walks by adopting a concept widely used in fluid 
dynamics [1], gas and plasma kinetics Q. For several 
models yielding normal and super-diffusive dynamics, we 
show that this function helps to uncover long-lived veloc- 
ity correlations that extend beyond the horizon dictated 
by the standard temporal correlation function. Further- 
more, the unfolding into the spatial domain allows for a 
meaningful description of velocity correlations even when 
the temporal correlation function docs not exist. We ar- 



gue that the process-specific generalized correlation func- 
tion can be accessed experimentally. Thus it can serve as 
a diagnostic tool to determine the underlying stochastic 
dynamics of real-life phenomena. 

Continuous-time random walks and s-t velocity auto- 
correlation function. We consider single-particle pro- 
cesses that belong to a class of continuous-time random 
walks (CTRWs) [31 • In its simplest one-dimensional re- 
alization, such a walk is performed by a particle moving 
ballistically with a fixed velocity Vi between two turning 
events. The duration of the j-th 'flight', that is the time 
interval between two consecutive turnings, = t^+i — t^, 
is governed by a probability density function (PDF) iP{t). 
At the end of each flight the particle changes its velocity 
to a new random value Wi+i sampled from the PDF h{v) 
and then starts a new flight. Two random variables, 
and Vi, are statistically independent but the spatial and 
temporal evolution of the walker during the flight is cou- 
pled, Xi+i —Xi = ViTi. This simple setup is able to repro- 
duce normal and anomalous diffusion by a proper choice 
of PDF's of flight-times and velocities [l^, There is 
a multitude of real-life systems and processes whose evo- 
lution in space and time can be described by this model. 
Examples range from the dynamics of tracers in turbu- 
lent flows [l^l and cold atoms in optical potentials [isl 
to foraging animals [3, [IB] , and migrating T-cells jl6j|. 

The key property of such CTRW models is a well- 
defined velocity of a walker at any instant of time. This 
feature allows us to introduce space-time velocity auto- 
correlation function for a walker's velocity by redefining 
the conventional expression d, Q , 



C{x,t) = {v{Q,0)v{x,t)) , 



(1) 



for a single-particle process. Namely, a particle starts its 
random walk with velocity vq = v{x = 0,t = 0). After 
time t the particle is found at a point x with some velocity 
v{x,t). Then the product of two velocities is averaged 
over all possible realizations of the random walk process 
and all possible initial velocities, v^, sampled with h{vo). 
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The definition can be rewritten in integral form, 



C{x,t) = / / vvoP{v,x,t\vo,0,0)h{vo)dvodv, (2) 



where P(u, x, t|wo, 0, 0) is the probability density func- 
tion of finding a particle in {x,t) with velocity v, pro- 
vided it has started with velocity vq at time t = and 
X = 0. Note that the function given by Eqs. (|ll2p . 
when being integrated over x, yields the standard ve- 
locity auto-correlation function C{t) ~ {v{0)v{t)) [l7| . 
In full analogy with the common derivation of the trans- 
port equations for the density of particles [loj . we can 
derive an equation for x, t|uo, 0, 0). We introduce 
the frequency of velocity changes, i'vg{x^t), which counts 
the number of particles whose flights ended in the in- 
terval [x, X + dx] during the time interval [t, t + dt]. The 
additional subscript vq tracks the history of particles and 
denotes only those which had velocity at t = 0. The 
balance equation for assumes the form of an integral 
equation, reading 

+ 00 t 

Vy^ix^t) ~ j '^^ j ^Vo{x — VT^t — T)h{v)ll){T)dT 
-oo 

+ ^p{t)S{x -vot). (3) 
A particle changes its velocity at the end of the fiight of 



duration r that was initiated at the point x — vt. Mul- 
tiplication by 'ip{T)dT and h{v)dv yields the probability 
of having a flight-time r and velocity v. We also assume 
that all particles start their random walk at t = and 
X = 0. The last term on the right hand side accounts for 
the particles that finish their very first fiight at the given 
moment of time t. Correspondingly, P(u, x, t|vo, 0, 0) is: 

t 

P{v,X,t\vo,0,0) ^ J Vygix - VT,t - T)h{v)^{T)dT 


+ 'f{t)5{x~vomv~vo). (4) 

A particle has a velocity v at point (x, t) if it has previ- 
ously changed its velocity at time t — t and still is in the 
process of flight with velocity v. The probability to stay 
in the flight until time t is given hy '^{t) = 1 — J* %l){T)dT. 
The second term on the right hand side of Eq. ^ ac- 
counts for the particles that are still in their first flight. 
As before, we track only the particles that started out 
with velocity vq. The above two equations can be re- 
solved by using a combined Fourier-Laplace transform 
with respect to x and t which turns the convolution-type 
integrals into algebraic products. We use C[..] to denote 
the Laplace transform, a hat for Fourier transforms, and 
a tilde for a combination of the two, whereas k and s de- 
note coordinates in the Fourier and Laplace spaces. We 
then flnd P{v,vq, k, s), and by using Eq.([2]), obtain the 
general expression for the s-t auto-correlation function in 
the Fourier-Laplace space. 
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(5) 



Equation ^ constitutes the central result of this Let- 
ter. An important observation can be made directly from 
this general expression: The integral of C{x,t) with re- 
spect to X, that corresponds to setting fc = 0, yields 
the standard velocity auto-correlation function C{t). Re- 
markably, the flrst term in the equation that describes the 
contribution from particles that changed their velocities 
several times, integrates to zero. It is only the particles 
remaining in their flrst flight during the time t that con- 
tribute to the C{t). This implies that C{t) (wg) 

Next we consider several particular models of the 
CTRW, covering the broad spectrum of diffusive regimes, 
from the standard diffusion to super-ballistic super- 
diffusion. 

Levy walks: from normal to ballistic super- diffusion. 
The Levy walk process is a particular case with 
a bimodal choice of the velocity PDF, h{v) — 
[d{v - vo) + d{v + vq)] /2, and the PDF of flight-time of 



a form 



(6) 



T-o(l + r/ro)i+^ ' 

where tq sets the time-scale of the process, and the scal- 
ing exponent is positive, 7 > 0. The latter parameter 
plays a key role in deflning the type of the diffusion. If 
7 > 2 then the mean square of the flight-time is finite 
and the process reproduces normal, Brownian-like dif- 
fusion, with the linear scaling of the mean squared dis- 
placement, {x^{t)) oc t 0. The mean squared flight-time 
diverges for 1 < 7 < 2, which leads to anomalously fast 
spreading, {x^{t)) oc t^~^ 0, [I!]- In the extreme case of 
< 7 < 1, the anomaly is strong and prevails the evo- 
lution, because even the flrst moment of the flight-time 
diverges; so that the mean squared displacement scales 
ballistically, {x'^it)) (x t'^. In all cases, however, the den- 
sity of particles P{x, t) is confined to the ballistic cone, 
X G [-vot,vot]. 
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When 7 > 2, in the inside of the casual cone and in 
the asymptotic hmit, the density obeys the standard dif- 
fusion equation Q , 



dP{x,t) 
dt 



= DAP{x,t); D 



7-2 



(7) 



In this case, the first term in the general formula ([5|) 
reduces to the simple result, denoted by 



C'centr(2;, t) — V^^ 



Dtp 

7-1 



AP{x,t) 



, To dPix,t) 

'7-1 dt ■ 



(8) 



Equation (|8]) reveals an interesting relation between the 
particle density function and the corresponding s-t cor- 
relation function, namely that C{x,t) is proportional to 
dP{x,t)/dt. Note that the contribution of the second 
term in Eq. ([5]) corresponds to the ballistic delta-peaks, 
running with the speed vq and decaying in time accord- 
ing to ^'(t). Ballistic peaks arc the hallmark of Levy 
walks [1^. Their contribution is typically considered to 
be asymptotically vanishing in the regime of standard dif- 
fusion. However, only the peaks contribute to the tempo- 
ral velocity correlation function C{t) and therefore they 
cannot be neglected. The results presented by Eqs. ([71 [S]) 
arc valid for an arbitrary choice of 'ip{t) that has a finite 
second moment, T'^'4){T)dT < 00, including the most 
generic form, an exponential function iP{t) = e~'^ 1'^° Itq. 
In this case C(t) = WqC"*/'^" while the space-time veloc- 
ity correlations, for example, at the point a; = scale like 
^-3/2^ When the flight-time is constant, ^/'(r) = (5(r — tq), 
C(t) drops to zero outside the region t < tq. At the same 
time, by keeping the spatial-dependence in the correla- 
tion function, Eq. ((!]), one can observe temporal corre- 
lations which extend outside this time-window. These 
two examples highlight the fact that s-t velocity corre- 
lation function provides access to long-lived correlations 
and therefore increases the chance of their detection. 

For 1 < 7 < 2 the mean squared flight-time diverges. 
It induces a super-diffusive behavior with the den sity of 
particles obeying a generalized diffusion equation [20| , 



dP{x,t) 
dt 



-K{~Ar/^P{x,t); 



where K = t^~\2 il - l)r[l - 7] cos(7r7/2) and A^/^ 
is the fractional Laplacian operator [2l|. Note that this 
description is valid in the inner part of the casual cone 
only. In there C{x,t) is proportional to the fractional 
Laplacian of the density of particles. 



C'contr(2;5 t) 



-vlIl^(-Ay/^P{x,t) 

7-1 
w^To dP{x,t) 
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1 



dt 



(10) 



Again, by virtue of the Eq. ([9]), in the asymptotic limit 
the central part of the s-t correlation function is given by 
the time derivative of the corresponding density. There- 
fore, the velocity auto-corrclations are negative near the 




FIG. 1. (color online) Space-time velocity auto-correlation 
function C{x,t) for the super-diffusive Levy walk, Eq. ((6)1, 
with 7 = 3/2 as a function of x and t. The red dashed lines 
indicate the positions of local maxima xjj^ on the x — t plane 
which follow the power-law scaling oc ±1^^'' , while the 
height of the maxima decays as t~^~^^'' . The inset depicts 
spatial profiles of C{x, t) at two different instants of time, 
t = 20 (heavy blue line) and 50 (light red line). 



point a; = 0, see in Fig. 1. Upon departure from the 
origin the correlations become positive and produce two 
local maxima. As before, the integral of the central part 
of the s-t function with respect to x vanishes, and only 
ballistic peaks contribute to G{t) in this case alike. 

By setting 7 < 1 in Eq. ©, one can enhance the 
anomalous character of the process. The average flight- 
time of particles diverges and this implies a ballistic scal- 
ing for the density of particles. The correlation C{x,t) 
can be evaluated in the Fourier-Laplace space. As an 
illustration, we consider the case 7 = 1/2 where both 
quantities, P{x, t) and C(x, tV can be expressed in terms 
of analytic functions [l3, [lH- The density exhibits an 
tZ-shaped proflle, diverging at the ballistic fronts: 



Pix,t) 



0{VQt 



nit^vl - a;2)i/2 



(11) 



The s-t velocity auto-correlation function behaves simi- 
larly: 



(9) C{x,t) 



2 1/2 



S{x — vot) S{x + vot) 



(12) 

eivot-\x\) 



<l/2 



^1/2 



2i3/2t>o 



From this result it follows that the s-t velocity correla- 
tion function is negative and nearly constant inside the 
ballistic cone. The relative decay rate P{x,t)/C{x,t) is 
of the order t^/^ now. The profiles evaluated via direct 
numerical simulations of the random walk, cf. Fig. 2(a), 
perfectly match the analytical prediction. 

Super-ballistic super- diffusion. Regimes of diffusion 
described by the Levy walk model are bound from above 
by the ballistic propagation so that no particles can cross 
the front \x\ = vot. One possible way to overcome this 
limitation is to allow the flight speed to have a broad 
distribution. Assume that h{v) is a Cauchy (also known 
as Lorentzian) distribution, h{v) = Vq^{1 -\- v'^/vq)~^, a 
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FIG. 2. (color online) Densities P{x, t) (dashed red lines) and 
space-time velocity auto-correlation functions C{x,t) (solid 
blue lines) of (a) the Levy walk in the regime of ballistic 
super-diffusion, Eq. with 7 = 1/2 and (b) for the super- 
ballistic super-diffusion with tp{T) = [S{t — 1) -f S{t + l)]/2 
(see text). Correlations functions were obtained numerically, 
by averaging over 10* realizations. The inset shows C(0, t) 
at four instants of time. The solid line corresponds to the 
power-law t^^^^. All other parameters are set to one. 



velocity PDF frequently employed in plasma and kinetic 
theories [l^. For such velocity PDF, the density of par- 
ticles is independent of the flight-time distribution and 
also possesses the Lorentzian shape in the asymptotic 
limit [13, 



vot 



(13) 



The expression for the s-t correlation function then ac- 
quires the form: 



Vot 



1 



7r(wQt2 + x"^) Vot 



(14) 



The first term is proportional to the density but with the 
opposite sign. The second term depends on the flight- 
time distribution and vanishes in the asymptotic limit. 
Note that the integral of C{x, t) vifith respect to the coor- 
dinate diverges meaning an infinite time-dependent cor- 
relation function C{t). In clear contrast, the s-t correla- 
tion function is well-defined, see Fig. 2(b). 

Discussion. A key feature of all considered random 
walks is the region of negative correlations near the vicin- 
ity of the starting point. This means that the major- 
ity of particles found there are flying backwards, in the 
direction opposite to the one of their original motion. 
This quantitative result may appear rather obvious when 
stated here; - yet, to the best of our knowledge, it has 
neither been presented nor proven before. The shape of 
the negative-correlation region and the time-scaling of its 
length and depth at a; = 0, are model-speciflc [2j]. Al- 
though we considered only a few models, the presence 
of strong negative velocity correlations around the point 
X ~ should hold true for any generic random walk 
process with flnite velocities of diffusing particles. We 
focused on the situation where all particles started their 



walks at time t = 0, but there is no conceptual difficulty 
to consider particles with an existing pre-history, or to 
address the case of initial velocity distribution different 
from the equilibrium one h(v). Moreover, velocity fluc- 
tuations during single flight events can also be included 
into the model 
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Perhaps the best candidate for the analog simulation 
of super-diffusive continuous-time random walks is a cold 
atom moving in a periodic optical potential. There are 
strong evidences, both theoretical [l^ [2^ and experi- 
mental [2^ [131 , that the dynamics of the atom can be re- 
produced either with the Levy-walk model or some mod- 
iflcations of it [13, ■ There exist several possibilities to 
prepare a strongly localized ensemble with all atoms hav- 
ing near equal velocities. They include either the sudden 
release of atoms from a ballisticaly moving deep optical 
well or more exotic setups like atomic lasers [29| or the 
controllable soliton emission [s^l as in the case of dense, 
ultra-cold ensembles. For such initial states the calcu- 
lation of s-t velocity auto-correlation function C{x,t) is 
equivalent to flnding the PDF of the spatial velocity dis- 
tribution P{v, X, t), see Eq. The PDF of velocities in- 
tegrated over spatial coordinate, P(w, t), can be measured 
by using the routine time-of- flight technique [3l|, when 
the instantaneous velocity of a released atom is trans- 
formed into the atom position and the latter is measured 
by using the florescence effect. It is also feasible to mea- 
sure the instantaneous spatial distribution of atoms by 
using the florescence image of the cloud, however, with- 
out resolving cloud's velocity distribution [3l|. A mea- 
surement of the space-dependent PDF of velocities will 
require the use of the time-of-flight technique combined 
with a consecutive de-convolution of the fluorescence im- 
age [s^l ■ By knowing the spatial distribution of the cloud 
at time t (obtained from another experiment under the 
same conditions), the de-convolution transform can be 
performed numerically thus finally yielding C{x,t). 

Conclusions. The spatial dependence of the s-t spa- 
tiotemporal velocity auto-correlation function, Eq. ([l}, 
can be decomposed into two contributions. A first one 
is created by the particles which undergo several fiights 
before the observation time. The second originates from 
the particles that are still in their first fiiglit. For any 
random walk with finite average fiight-times, the central 
part of the velocity correlation pattern can be calculated 
as the time derivative of the particle's density P{x,t). 
Noteworthy, the generalized diffusion equations for the 
density of particles, Eqs. (|7l9p . do not contain the veloc- 
ity degree as a dynamical variable; only its absolute value 
enters the equations indirectly, via the diffusion constant. 

We believe that the concept developed here can be uti- 
lized for any process that is described as a continuous- 
time random walk and where a physically finite veloci- 
ties can be attributed to the diffusing entities at any in- 
stant of time (33l - [35| . The spatiotemporal velocity auto- 
correlation function can be generalized to the case of 
two-dimensional random walks and thus provides another 
level-of-detail for the analysis of the transport phenom- 
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ena from different areas of physics and real life that can [ly, [2^, [Sa] . 

be modeled with continuous time random walks [l^, [13- This work has been supported by the DFG Grant 

HA1517/31-2 (S.D. and P.H.). 
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